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Abstract 

Under the combined constraints of rapid rotation, sphericity, and magnetic field, 
motions in planetary cores get organized in a peculiar way. Classical hydrodynamic 
turbulence is not present, but turbulent motions can take place under the action 
of the buoyancy and Lorentz forces. Laboratory experiments, such as the rotating 
spherical magnetic Couette Z? 75* experiment in Grenoble, help us understand what 
motions take place in planetary core conditions. To cite this article: H.-G. Nataf 
and N. Gagniere, G. R. Physique XXX (2008). 

Resume 

Sur la nature particuliere de la turbulence dans les noyaux planetaires. 

Sous les contraintes combinees de la rotation rapide, de la sphericite et du champ 
magnetique, les ecoulements dans les noyaux planetaires s'organisent d'une maniere 
particuliere. La turbulence hydrodynamique classique n'est pas presente mais des 
mouvements turbulents peuvent se mettre en place sous Paction des forces d'Ar- 
chimede et de Lorentz. Des experiences de laboratoire comme I'experience DTS de 
Couette spherique sous champ magnetique a Grenoble, nous aide a comprendre les 
ecoulements qui peuvent exister dans les conditions des noyaux planetaires. Pour 
citer cet article : H.-G. Nataf and N. Gagniere, G. R. Physique XXX (2008). 
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1 Introduction 



The recent success of the VKS dynamo experiment [1,2], after the harvest of 
the pioneer experiments in Riga [3] and Karlsruhe [4], brings hope to better 
understand the mechanisms at work in the dynamo process. It is believed that 
the self-sustained dynamo process generates the magnetic fields of most as- 
trophysical objects, including our planet, the Earth. In this process, a given 
velocity field u can produce a magnetic field controlled by the magnetic in- 
duction equation (1), which shows that if the velocity is large enough for the 
induction term V x (jl x to dominate over the diffusion term rjAB, the 

magnetic field B can grow. 

= ^x (uxB) +r]AB (1) 
ot ^ ^ 

The magnetic field of the Earth originates in its core, made of liquid iron. 
As the Earth cools down, on a geological time-scale, heat is carried from the 
core to the mantle by convective motions. These motions are governed by 
the Navier-Stokes equation, which, in the Boussinesq-approximation may be 
written as: 

fdu \ - - ^ ^ ^ 

p ( — + n- Vnj = -VP - 2pn x u + fiAu + j x B + fa (2) 

where the various symbols have their usual meaning. On the right hand side, 
one recognizes the Coriolis force —2Q x u, the viscous force nAu, the Lorentz 
force j X B, and the buoyancy force fa- 

When plugging in typical values of the properties of liquid iron and of the ve- 
locity u and magnetic B fields, one finds that, because of the fast rotation Q, 
the dominant forces are the Coriolis and Lorentz forces. The resulting regime 
is called "magnetostrophic" , and it has received the attention of theoreticians 
for decades. Taylor [5] was the first to point out that, in this regime, if one 
considers tubes of liquid co-axial with the rotation axis, there is no torque 
to resist the torque applied by the Lorentz forces. He concluded that, in the 
steady-state, the torque applied by the Lorentz forces must vanish, yielding 
what is now called a "Taylor-state". The magnetostrophic regime and the 
expected Taylor-state remain very difficult to reach in numerical simulations 
(because it is technically difficult to neglect viscous terms) and their explo- 
ration in laboratory experiments is recent [6,7,8]. 

In this article, we will highlight some striking properties of the "magne- 
tostrophic" regime, as derived from measurements performed in the DTS" ex- 
periment. In section 2, we describe the DTS setup and present the relevant 
parameters. The time-averaged flow is analyzed in section 3, and the (weak) 
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Figure 1. Sketch of the central part of the DTS experiment. The 7.4cm-radius inner 
sphere is made of copper and contains 5 layers of magnets. The inner radius of the 
stainless steel outer sphere is 21cm and its thickness is 5mm. Blind holes are drilled 
and threaded along several meridians and parallels to receive electrodes. The large 
holes at latitudes —20°, 10° and ±40° receive interchangeable assemblies, which can 
be equipped with ultrasonic transducers. 

fluctuations in section 4. Implications for turbulence in planetary cores is dis- 
cussed in section 5. 



2 The DTS experiment: exploring the magnetostrophic regime 

The Derviche Tourneur Sodium (DTS) experiment has been designed for the 
exploration of the magnetostrophic regime, in which the Coriolis and Lorentz 
forces are dominant [9]. It consists in a spherical Couette flow where both 
the inner and outer spheres can rotate rapidly at separate angular velocities 
around a common vertical axis. Forty liters of liquid sodium fill the shell 
between the two spheres and the inner sphere contains a strong magnet. The 
set-up is depicted in figure 1. The fluid flow is thus strongly influenced by both 
the Coriolis force and the Lorentz force. Dimensions and typical dimensionless 
numbers are given in Table 1. 

The radius of the outer sphere is a = 21cm. Because of this moderate di- 
mension, the Joule dissipation time tj of magnetic fluctuations by diffusion is 
small (less than one tenth of a second). The Alfven wave velocity Ua depends 
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Table 1 

Typical values of the relevant parameters and dimensionless numbers for given im- 
posed rotation frequencies / = Q/27r of the outer sphere and differential rotation 
of the inner sphere A/ = AO/27r with respect to the outer sphere. For the num- 
bers that depend on the magnetic field strength, two values are given, the first one 
with B = Bi = 0.175 T at the equator of the inner sphere, the second one with 
B = Bo = 0.008 T at the equator of the outer sphere. 

on the strength of the magnetic field, and is thus much stronger near the 
inner sphere than at the outer sphere. Comparing the Alfven time to the dif- 
fusion time yields the Lundquist number S. It reaches 12 on the inner sphere, 
suggesting that Alfven waves can survive there, while they must be severely 
damped near the outer boundary, even though the Hartmann number Ha is 
large everywhere in the shell. 

The rotation frequency / of the outer sphere has been varied between and 
15Hz. The small value of the Ekman number for a typical rotation frequency 
of 5Hz illustrates that the Coriolis force dominates over viscous forces in the 
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bulk of the fluid. A very thin Ekman layer (one tenth of a millimeter) must 
form beneath the outer shell. The value of the Elsasser number A, which 
compares the Lorentz force to the Coriolis force, shows that rotation effects 
dominate in that region, while the opposite holds near the inner sphere. The 
domain of magnetostrophic equilibrium A ~ 1 occupies the central part of 
the fluid shell. The number A, called Lehnert number in [10], compares the 
frequencies of Alfven modes to that of rotation (inertial) modes. It is smaller 
than 1 everywhere. 

Other effects depend upon the strength of the forcing, measured by A/, the 
differential rotation frequency of the inner sphere with respect to the outer 
sphere. It can be varied between about — 20Hz and 20Hz. For a moderate value 
of 5Hz, the expected fluid velocity U is of the order of 2m/s. The Reynolds 
number Re is thus very large, while the magnetic Reynolds number Rm is 
larger than 5, indicating that the magnetic fleld is modifled by the flow. The 
large value of the interaction parameter N near the inner sphere demonstrates 
that, in turn, the magnetic fleld deeply influences the flow. 

The DTS set-up is well suited for exploring the magnetostrophic regime. It is 
not a dynamo experiment, but the imposed magnetic fleld is strong enough to 
give rise to Lorentz forces of the same order as the Coriolis force. Both forces 
have a strong influence on the flow. The magnetic Reynolds number is large 
enough for induction effects to be well developed. We now examine what flow 
results from this complex combination, which is typical in planetary cores. 



3 Characteristics of the mean axisymmetric flow 

We use ultrasonic Doppler velocimetry to measure the mean axisymmetric 
flow. A transducer is placed in one of the interchangeable assemblies at latitude 
10°, and shoots a beam with declination and inclination angles of 60° and 
66°, respectively. After a 38cm-long linear trajectory, the beam hits the outer 
sphere at a latitude of about —28°. Azimuthal velocities being more than 10 
times larger than radial velocities, the along-beam measured velocity directly 
yields an angular velocity proflle. 

Such a proflle is shown in flgure 2 for / ~ 5Hz and A/ ~ 5Hz. Starting from 
the beam entry point at latitude 10°, we see that the angular velocity rises 
regularly from a value of about 0.25m/s, reaching an almost uniform plateau 
farther in the fluid, and increasing again to reach a maximum of 1.05m/s 
where the beam gets closest to the inner sphere. Comparing with axisymmetric 
numerical simulations, Nataf et al. [8] explain that the velocity bump near the 
inner sphere is a signature of the magnetic wind: the Elsasser number is large 
and the angular velocity isolines follow the magnetic fleld lines; the magnetic 
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coupling between the rotating inner sphere and the hquid sodium is efficient. 
Away from that bump, the flow is geostrophic as the Elsasser number gets 
below 1 (see Table 1). The geostrophic flow is driven by the magnetic torque 
that results from the shearing of the magnetic fleld lines and slowed down by 
friction in the Ekman layers beneath the outer sphere. This situation has been 
well studied by Kleeorin et al. [11], in the asymptotic limit of small Ekman, 
Rossby, Reynolds and Elsasser numbers. At small cylindrical radius s, the 
imposed magnetic fleld is large and friction in the Ekman layers of the outer 
sphere is weak. As a consequence, the shear must be very small in order for 
the magnetic torque to remain small enough to balance the friction torque: 
the fluid is in rigid body rotation. The regular decrease in angular velocity as 
one gets farther away from the inner sphere was predicted by Kleeorin et al. 
[11]: friction increases while the imposed magnetic fleld decreases, implying 
an increasing shear. It is important to realize that the flow is geostrophic, as 
a consequence of strong rotation, but entirely driven by magnetic forces. A 
modifled Taylor state is achieved, where the magnetic torque on cylindrical 
tubes parallel to the axis of rotation only balances the (weak) friction in 
the Ekman boundary layers. Note that, under this mechanism, the induced 
toroidal magnetic fleld is expected to be larger in the external region, where 
shear is large, than in the region near the inner sphere, although the imposed 
magnetic fleld is much larger there. 

We flnd similar proflles for all values of A/ as long as the Rossby number 
Ro = A/ / / remains smaller than a few units. A good quantitative flt between 
the theory and the measurements is achieved, when one takes into account that 
the Ekman layers become turbulent in the experimental conditions (see [8]). 
Therefore, we think that the same dynamics prevails for the very small Rossby 
numbers that characterize flow in the Earth's core. 



4 Weak turbulence 

The proflle shown in flgure 2 is, in fact, a histogram constructed from over 1000 
shots. Although the Reynolds number reaches 7- 10^ for this run (see Table 1), 
the fluctuations around the average proflle are amazingly small. Besides, the 
power spectra of the fluctuations reveal peculiar bumps at several frequencies. 
This is best seen on the spectrogram of a whole run, where / ^ 5Hz and A/ 
is varied in steps from to — 20Hz. 

Here, the signal is the difference in electric potential between two electrodes 
at the same latitude of —45° and distant of 40° in longitude (see flgure 1). The 
spectrogram is shown in flgure 3. Frequency is on the horizontal axis and time 
on the vertical axis. The power spectral density of the signal is color-coded 
(log-scale). Thin lines reflect the two forcing frequencies / and — A/ and some 
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Figure 2. Profile of tlie angular velocity (velocity along the beam) of the flow as 
a function of distance from the outer sphere, obtained by ultrasonic Doppler ve- 
locimetry. The angular velocity rises smoothly to a plateau before it increases again 
as the ultrasonic beam gets closest to the inner sphere. / ~ 5Hz and A/ ~ 5Hz. 
The (non-linear) scale on the right gives the value of the adimensional cylindrical 
radius s. 

overtones. For —A/ < /, the spectrum is featureless, while bands of higher 
power are clearly visible when A/ reaches — lOHz. The detailed study of these 
bands performed by Schmitt et al [7] reveal that they correspond to waves 
that propagate in the direction of the flow (in the frame of reference rotating 
with the outer sphere) but at a slower velocity. The lowest band corresponds 
to a wave or mode with an azimuthal wave-number m = 1, the following 
one to m = 2, and so on. The origin of these waves is not clear yet. They are 
different from the inertial modes identified by Kelley et al [12] in a similar set- 
up but with a weak magnetic field. Indeed, the derivation of their dispersion 
relationship by Schmitt et al [7] reveal that their frequencies can be higher 
than that of inertial modes (which are bounded by 2 f fluid, where f fluid is the 
rotation frequency of the fluid in the laboratory frame of reference). 

Once again, useful information can be gained from ultrasonic Doppler velocity 
measurements. The bands that were first seen on power spectra of the induced 
magnetic field are also present in the power spectra of velocity, as illustrated 
in figure 4. They clearly extend deep into the sphere. 

We have seen that, in the modified Taylor state, the maximum induced toroidal 
field is not near the inner sphere. However, it is clear that a local perturbation 
of the azimuthal velocity near the inner sphere will induce a strong azimuthal 
magnetic field, which will in turn change the magnetic field and the torque 
that drives the azimuthal flow. It is tempting to see in this effect the origin of 
the waves we observe in DTS, but further work is needed to corroborate this 
hypothesis. 
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Figure 3. Spectrogram of electric potential measurement. The signal is the electric 
potential measured between two electrodes situated at latitude —45° and distant 
of 40° in longitude. The colors give the log-amplitude of the short-time spectral 
density. The horizontal axis is frequency and the vertical axis is time. During this 
time-frame, the outer sphere was rotating at a constant angular frequency / = 5Hz, 
while the inner sphere differential rotation rate A/ was varied in steps from to 
— 20Hz. The thin lines, labelled / and A/, reveal these frequencies and some over- 
tones. The spectrum is almost void until A/ reaches about — lOHz. The spectrum is 
then dominated by 3 broad bands (Wl, W2, W3) whose frequency slightly decreases 
as I A/ 1 increases. 
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Figure 4. Radial profile of power spectral density of radial velocity versus frequency. 
Radial velocities are measured by ultrasonic Doppler velocimetry along a radial 
beam shot from a transducer at —40° in latitude. We then derive this image of 
color-coded log of power-spectral density as a function of adimensional radius r/a 
and frequency. The rotation frequency of the sphere is / = 4.3Hz and the differential 
rotation frequency of the inner sphere is A/ = — 23Hz. The power is high in bands 
at several frequencies, as in the spectrogram of figure 3. 

Radial velocity spectrum 
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Figure 5. Power spectral density of radial velocity versus frequency at mid-depth. 
For the same record as in figure 4, we plot the log of power spectrum measured at 
mid-depth. Peaks at several frequencies dominate the spectrum. 
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5 Implications for turbulence in the Earth's core 

Our experiments illustrate that, due to the strong influence of the Coriolis 
force, motions are geostrophic when the Elsasser number is less than 1, even 
though the magnetic torque on co-axial tubes is the only driving mechanism. 
Jault [10] argues that this property remains true even for large Elsasser num- 
bers if one considers time scales that are short compared to the Alfven time. 
Pais and Jault [13] apply this idea to retrieve geostrophic motions in the core 
that explain the observed secular variation of the magnetic fleld of the Earth 
between years 2000 and 2005. 

The idea that motions are organized in cylinders around the axis of rotation 
dates back to Taylor [5]. He pointed out that, in the steady-state, the mag- 
netic torque acting on such cylinders must vanish since, in the limit of small 
Ekman and Rossby numbers, there is no torque to balance it. The DTS ex- 
periment supports this idea, with the modiflcation that the magnetic torque 
is balanced by friction at the core-mantle boundary [11], and that this fric- 
tion takes into account that the Ekman boundary layer becomes turbulent 
[8]. Although the local Reynolds number of the Ekman layers remain small in 
the core [14], enhanced friction at the core-mantle boundary can result from 
a rough boundary [15,16], or from electromagnetic coupling between the core 
and the mantle [16]. In any case, nutation observations require a 10000-fold 
enhancement of core-mantle coupling, as compared to linear Ekman friction 
[17]. It would be worth it to take this effect into account in numerical simu- 
lations, in particular in quasi-geostrophic dynamo simulations [18], which are 
particularly efficient for exploring Earth-like fluid properties, and in which 
Ekman friction is already parameterized. 

Our experiments also suggest that the only fluctuations left in the magne- 
tostrophic regime, apart from localized turbulence in the Ekman layers, are 
waves that propagate azimuthally in a retrograde fashion. These global wave 
phenomena are usually excluded from the local turbulence analyses used to 
calibrate sub-grid algorithms in 3D numerical simulations [19]. This should 
probably be re-assessed once the origin of these waves is clearly established. 
In that respect, one should note one important limitation of our DTS experi- 
ment: the time-scales that characterize the various kinds of waves are all very 
similar. For typical values of / and A/, the period of inertial waves and Alfven 
waves, as well as the dissipation time of Alfven waves are all of the order of 
0.1s. This makes it difficult to clearly identify the origin of the waves, and it 
prevents a complete exploration of the various regimes. 

It was soon recognized that departures from the Taylor state would trigger 
a special kind of Alfven waves, called torsional oscillations [20]. Assuming 
reasonable values for the -hidden- magnetic fleld within the core, the period 
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of these waves should be of the order of 10 years. Both core flow models derived 
from the secular variation of the magnetic field [21] and numerical simulations 
of the geodynamo [22] have been explored in order to find evidence for these 
oscillations. The results indicate that part of the fluid motions could indeed 
be torsional oscillations, providing a nice link with the observed length-of-day 
variations [23] . However, recent analyses of the more precise satellite data also 
reveal large-scale motions that are not torsional oscillations. Regional "jerks" 
are thus observed by Olsen and Mandea [24], while the quasi-geostrophic core- 
flow models of Pais and Jault [13] are characterized by a large retrograde jet 
that is completely offset with respect to the axis of rotation, in violation of 
the axisymmetric character of torsional oscillations. All these observations, 
together with our detection of non-axisymmetric propagating waves in the 
DTS experiment suggest to us that the importance of torsional oscillations 
in the Earth's core could have been overemphasized. However, we should note 
again that, because the Alfven time in our experiment is of the same order as 
the Joule dissipation time, torsional oscillations must be severely damped. It 
would be of great value to conduct experiments where Joule dissipation is less 
important. 

Finally, we emphasize that fluctuations are very minute when the flow is un- 
der the combined influence of a strong magnetic field and rapid rotation (see 
figure 2 and 3). Basically, the flow is entirely constrained once the magnetic 
field and the rotation are given. When the Elsasser number is small, or when 
the timescale of the motions is short [10], the motions are constrained to be 
invariant along the rotation axis (quasi-geostrophy) but are driven by the 
Lorentz forces. As magnetic diffusion is small at short time-scales, only the 
flow can modify the magnetic field, and the same holds for the buoyancy field. 
One is thus led to conclude that the turbulence at work in planetary cores has 
little to do with classical turbulence. Our DTS experiment is missing some 
of the key ingredients of a planetary dynamo, in fact it is not a dynamo. The 
VKS experiment, which does behave as a dynamo [1], is also missing key in- 
gredients of planetary dynamos: turbulence is almost entirely hydrodynamic 
and the produced magnetic field has a negligible effect on the flow. The mech- 
anisms of the spectacular magnetic reversals observed in VKS [2], as those of 
the pioneer experiments of Lowes and Wilkinson [25] , probably have little to 
do with those at work in the Earth's core. 



6 Conclusion 

Laboratory experiments such as our DTS experiment help us better under- 
stand the peculiar organization of fluid flow and magnetic field at work in 
planetary dynamos. When the Coriolis and Lorentz forces are dominant, and 
the Elsasser number less than 1, a modified Taylor state is observed, where 
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the geometry is given by the rotation (invariance along the axis of rotation) 
and the driving by the Lorentz force. In a Taylor state, the magnetic torque on 
co-axial tubes self-adjusts to zero, since there is no balancing torque. In our 
experiments, as foreseen by Kleeorin et al [11], Ekman friction at the external 
surface balances the magnetic torque. In this modified Taylor state, the flow 
remains geostrophic (azimuthal) and its variation with cylindrical radius is 
entirely controlled by the balance between the magnetic torque the resulting 
shear produces and the friction at the external boundaries. Using ultrasonic 
Doppler velocimetry, we measure time-averaged azimuthal velocity proflles in 
perfect agreement with this analysis. 

Under the combined constraints of the imposed rotation and magnetic fleld, 
velocity fluctuations are very minute even though the Reynolds number is in 
the range 10^. When present, the fluctuations reveal a very peculiar behavior: 
the power spectra are dominated by peaks at several frequencies. The detailed 
analysis of the signals reveal that they correspond to waves that propagate in a 
retrograde fashion with respect to the fluid [7]. The spectral peaks correspond 
to various azimuthal wave numbers. It is not yet clear what are the driving 
forces of these magneto-inertial waves. 

We infer that turbulence at work in planetary dynamos is of a very special 
kind: the flow is a complete slave of rotation and magnetic forces, Ekman 
friction at the boundaries providing a balancing torque. But the magnetic 
fleld itself is advected by the flow, and the same holds for the buoyancy fleld. 

Two main limitations of the DTS experiment prevent it of approaching regimes 
more relevant for planetary situations. The flrst one is that it is not a dynamo 
experiment, hence the magnetic fleld is not free to be advected by the flow. 
However, in most experimental dynamos, the produced magnetic fleld is too 
weak to substantially alter the flow, in contrast with what occurs in DTS. 
The second limitation is that all relevant time-scales (Alfven, Joule dissipa- 
tion, rotation, Rossby times) are in the same range (of order 0.1s), in contrast 
with planetary situations. This makes it very difficult to isolate which are the 
main forces controlling the propagating waves we observe, and extrapolating 
to natural situations. Building an experiment exempt of these two limitations 
represents a tantalizing challenge, which has led Dan Lathrop and his group 
to build a promising 3m-diameter experiment in Maryland. 
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